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tained using very well-known ^He-'^He and T4,-T4, interaction potentials and several 
models for the ^He-TJ, interatomic potential. All the calculations have been per- 
formed with variational and diffusion Monte Carlo methods. It takes at least three 
atoms to form a mixed bound state. In particular, for small clusters the binding ener- 
gies are significantly affected by the precise form of the ^He-TJ, interatomic potential 
but the stability limits remain unchanged. The only exception is the '^He2TJ, trimer 
whose stability in the case of the weakest ^He-TJ, interaction potential is uncertain, 
while it seems stable for other potentials. The mixed trimer ^He(T4.)2, a candidate 
for Borromean state, is not bound. All other studied clusters are stable. Some of 
the weakest bound clusters can be classified as quantum halos, as a consequence of 
having high probability of being in a classically forbidden region. 
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I. INTRODUCTION 



There are many similarities between spin-polarized tritium (TJ,) and ^He atoms. Both are 
bosons, have a small mass and a weakly attractive interaction potential at large distances. 
Just like ^He, bulk TJ, remains liquid in the limit of zero temperature and zero pressure. 
In 1976, Stwaley and Nosanow suggested that TJ, should behave very much like liquid 
'^He and therefore constitute a second example of bosonic superfluid.- Many similarities 
between the two systems have been confirmed by recent quantum Monte Carlo^ results. For 
instance, bulk TJ, is a superfluid at zero temperature, with a condensate fraction rao=0. 129(3) 
at the equilibrium density. Miller examined the bulk ^He-TJ, system as an example of 
binary mixtures using a variational ansatz and found that it prefers to be completely phase- 
separated.- 

^He clusters have been extensively studied, both theoretically and experimentally, ex- 
hibiting nanoscopic superfluidity. Due to the similarities between ^He and TJ,, a similar 
behavior of spin-polarized tritium clusters is expected. So far, TJ, clusters have been inves- 
tigated only theoretically. Small clusters of TJ, were for the first time studied by Blume et 
alA. The smallest bound TJ, cluster is the trimer, showing characteristics of a Borromean 
or halo state.-!^ Its tiny binding energy of only around 4 mK has been later confirmed by 
other authors^i^. The extension of the study to larger clusters^ confirmed that TJ, clusters 
are much more weakly bound and diffuse than the ^He clusters with the same number of 
atoms. The stability of mixed clusters of T4-->^ with spin-polarized hydrogen (HJ,) and 
spin-polarized deuterium (DJ,), in particular small ones, has also been investigated. For the 
clusters with DJ, the stability limits depend on the number of DJ, atoms and the occupation 
of its nuclear spin states. On the other hand, due to the small mass of HJ,, it has been shown 
that even 60 TJ, atoms are not enough to bind one HJ, atom. 

Mixed ^He-TJ, clusters have not been studied yet. TJ, has almost the same mass as ^He, 
but it is a boson, so it will exhibit different physical properties. Because of the small mass of 
both ^He and TJ, and weak attractive parts of the interaction potentials it is expected that 
the smallest clusters could be extremely weakly bound. If ^He(T4,)2 is bound, it would be 
an example of yet another Borromean state, because none of its subsystems are bound. The 
only mixed molecular system for which Borromean state has been predicted theoretically so 
far is ^He2K.— Recently, small mixed helium- hydrogen clusters ^He2H, ^He2H~, '^HeH2 have 

2 



been studied as well.— 

Weakly bound clusters are also candidates for Efimov states,— which are predicted to 
exist also in mixed systems. Interest in these states and other universal binding properties 
of small clusters has been significantly intensified in the last couple of years by the research in 
the field of cold gases. After the experimental detection of an Efimov state in an ultracold gas 
of cesium atoms,— giant Efimov trimers have been detected in other systems^^i"— including 
those of mixed species.— Later work was devoted to the study of the Efimov spectrun>2^ as 
well as to the predictions^ and detection^^ of giant tetramers. Furthermore, a recent study 
of weakly bound bosonic clusters has found a series of universal cluster states, that can be 
qualitatively interpreted by adding one particle at a time to an Efimov trimer.— 

Although the physics of molecular weakly bound clusters is not the same as the one of 
ultracold gases, it is interesting to investigate if similar behavior of few-body states appear, 
indicating the universality. Universal physics is also connected to physics of quantum ha- 
los,— "— defined as bound states of clusters of particles with a radius extending well into the 
classically forbidden region. It is therefore worth studying which small molecular clusters 
qualify as quantum halos. 

In this work, we report the ground-state energy of clusters having up to 8 ^He and 4 
Tl atoms, as well as their structural properties, obtained with the diffusion Monte Carlo 
(DMC) method. In Sec. II, we report briefly the DMC method and discuss the trial wave 
functions used for importance sampling. Sec. Ill reports the results obtained by the DMC 
simulations. Finally, Sec. IV comprises a summary of the work and an account of the main 
conclusions. 

II. METHOD 

A. Interaction potentials 

We have modeled the ^He-^He interactions with the Aziz HFD-B(He) interaction poten- 
tial^. Some calculations have also been performed with the Korona et al. SAPT potential.— 

The interatomic interaction between tritium atoms is described with the spin-independent 
central triplet pair potential which was determined in an essentially exact way by Ko- 

los and Wolniewicz.— As in our recent DMC calculations of bulk HJ, and TJ,,-i^ we have 
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used the recent extension of Kolos and Wolniewicz data to larger interparticle distances by 
Jamieson et al. (JDW).— The potential is finally constructed using a cubic spline interpo- 
lation of JDW data, which is smoothly connected to the long-range behavior of the T4,-T4, 
potential as calculated by Yan et al— The JDW potential used in the present work has 
a core diameter a = 3.67 A and a minimum of —6.49 K at a distance 4.14 A. We have 
previously verified that the addition of mass- dependent adiabatic corrections (as calculated 
by Kolos and Rychlewski^^) to the JDW potential does not change the energy of bulk spin- 
polarized tritium.- It is worth mentioning that within the Born-Oppenheimer approximation 
it has been explicitly shown that in the spin-aligned electronic state, tritium nuclei behave 
as effective bosons.— 

Several forms of ^He-TJ, interaction potentials are available in the literature, including 
those resulting from ab initio calculations and semiempirical potentials. In 1984, Jochemsen 
et a/.— proposed a semiempirical potential (R2) after discussing available potentials until 
that time and comparing theoretical data with diffusion experiments at low temperature. 
That potential is essentially the same as an older potential by Das et alM (DWW), which 



the authors of Ref. |38| themselves preferred to use. The most sophisticated ab initio cal- 
culation is that of Meyer and Frommhold from 1994 (MF), who obtained good agreement 
with experiment, except at very low temperatures.— Subsequently, Chung and Dalgarnc4i 
introduced slight modifications in the short-range repulsive part of the MF potential that 
resulted in a better agreement with diffusion measurements at low temperatures (MFmod). 
In both MF and MFmod potentials, the construction of the potential-energy curve at large 
separations and the dispersion component that determines the long-range behavior in MF 
and MFmod is of high precision.— All of our calculations have been performed with both the 
DWW and the MFmod potential. In addition, for some selected clusters we have performed 
calculations using some other forms of available potentials. Specifically, Toennies et alM 
have proposed the Lennard- Jones 6-12 potential (TWW) whose parameters have been de- 
termined using data from low energy elastic scattering. TWW has the largest core, smallest 
depth and the strongest long-range part of all the studied potentials. Due to its simplicity it 
has been used also in some very recent calculations, e.g. by Krotscheck and Zillich;— these 
authors have found that, in the bulk liquid, the binding energies obtained with TWW differ 
less than 5% from those that are obtained with the DWW potential. Finally, we include 
the MF potential and the potential by Tang and Yang^ (TY), who used the Tang- Toennies 
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FIG. 1. Comparison among interaction potentials between '^He-^He (HFD-B(He>22.), ''He-T| (TY-M, 
TWW-^, DWW2^, MFmod^'^) and T|-Ti (JDW-M). For long distances, TY and MFmod have 
basically the same behavior. 

model. The comparison of these different potentials is presented in Figure [1] The strongest 
MFmod potential has a core diameter of cr = 3.10 A and a minimum of —7.14 K at a 
distance 3.52 A. The MF potential would be the same in the scale of the figure because 
the only differences are in the repulsive part. The DWW interaction has a core diameter 
of cr = 3.18 A and a minimum of —6.53 K at a distance 3.60 A, while the weakest TWW 
potential has a core diameter of a = 3.31 A and a minimum of —5.34 K at a distance 3.72 A. 
In Fig. [T]the ^He-^He and T4,-T4, interaction potentials have been included for comparison 
as well. Of all the interaction potentials, the T4,-T4, one has the largest core and the most 
attractive long-range part. 




B. Diffusion Monte Carlo method 

The ground-state properties of the clusters have been studied using the DMC method, 
whose starting point is the many-body Schrodinger equation written in imaginary time, 

-h^^^ = iH-E,)^iR,t), (1) 

where Ej. is a constant acting as a reference energy and R= (ri, . . . , r^r) collectively denotes 
particle positions. 
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The A^-particle Hamiltonian H for the cluster '*He„(T4,)m, n + m = is 

N ,2 ^ ^ n+m n+m 

^ = - E + E ^^^eHe(r.,) + E E ^^^eT(r.,) + ^ l^Trlr,,) , (2) 

i=l *,i=l i=l j=n+l i,j=n+l 

i<j i<j 

where VHeHeir)^ yneTir) and VTrir) are the interaction potentials between the different 
components of the mixture. 

In DMC, the Schrodinger equation ([1]) is solved stochastically by multiplying \l/(-R, t) with 
ipiR), a trial wave function used for importance sampling, and rewriting Eq. ([1]) in terms 
of the mixed distribution t) = {R, t)ilj (R) . Within the Monte Carlo framework, 

$(-R, t) is represented by a set of walkers Ri. In the limit t — oo, only the lowest energy 
eigenfunction, not orthogonal to ip{R)i survives and then the sampling of the ground state 
is effectively achieved. Apart from statistical uncertainties, the energy of a A^-body bosonic 
system is exactly calculated. 

In the present simulations, the trial wave function has been written as a product of 
Jastrow two-body correlation functions between all the pairs, 

n n+m n n+m 

MR)=X{fHeHe{n,) H Mn,)^ n /HeT(r.,), (3) 
i,j=l i,j=n+l 1=1 j=n+l 

i<j i<j 

where fueUe describes two-body correlations between helium atoms, /tt describes two-body 
correlations between spin-polarized tritium atoms and fueT accounts for the "^He-TJ, pairs. 
Different types of two-body correlation function /(r) have been tested. 



/i(r) = exp 



sr 



(4) 



r /«\ T "I 

/2(r) = exp (^-j - sirj /r , (5) 
/3(r) = exp[-6i exp(-62r) - S2r\ (6) 

where 6, s, a, 7, Si, 61, 62 and S2 are variational parameters. The optimization of the 
trial wave functions has been done for all clusters by means of the variational Monte Carlo 
(VMC) method. For clusters having less than 5 atoms it has been very difficult to obtain 
good quality VMC results so different forms have been tried out to ascertain that DMC 
gives the same value of the energy for different guiding wave functions. 
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For ^He-^He correlations, a set of calculations for all clusters has been done using fi{r) 
dl]), with the parameter 6=2.6 A and s from 0.08 to 0.13 A~^. For some smaller clusters, 
f2{^) dSD has also been used, with optimal parameters in the range a G [2.75,2.82] A, 
7 G [4.1,4.7], si G [0.001,0.015] A^^ For Ti-Ti two-body correlations, the form f^ir) © 
has been chosen for all the clusters, based on the previous experience with pure and mixed 
TJ, clusters.- The optimal parameters have been obtained in the range: bi from 82 to 94, 62 
from 1.23 to 1.24 A~^, and S2 from 0.04 to 0.06 A~^, where smaller values of s correspond 
to larger clusters. The only exceptions are the clusters ^He2(TJ,)2 and ^He(TJ,)3^4 for which 
better energies were obtained with bi = 200, 62 = 1-5 A~^ and S2 = 0.065 A~^. In addition, 
for some of the smallest clusters variational calculations have also been done with the form 
/i(r) (jlj), where optimal parameters have been found around a = 4.2 A, 7 G [3.15, 3.82] and 
si ^ 0.001 k-\ 

The optimization of the ^He-T^ correlations has been the most demanding: it has been 
done for several potentials and forms of trial wave functions. First, for the DWW potential, 
/i(r) has been used with b around 3 A and s ranging from 0.04 to 0.06 A~^. With /i(r) and 
clusters with more than 5 atoms the VMC energies approach the 80% of the DMC values. 
Therefore, for the smallest clusters, /2(^) has been tried as well, with resulting parameters 
a G [3.59,3.69] A, 7 G [3.08,3.20], si = 0.001 A'^ In conclusion, VMC energies reached 
more than 40% of the DMC ones in all cases except ^He2T4,, where negative energy could 
not be obtained at the VMC level for the DWW potential. For this cluster, the best set of 
parameters for the function f2{r) is a = 3.62 A, 7 = 3.1, Si = 0.001 A~^, although it arrives 
to the same DMC value for other sets of parameters. In addition, ^He-TJ, correlations for 
that cluster have also been modeled with fsir) ([6]) (61 = 300, 62 = 1-93 A~^ and S2 = 0.05 
A~^) and despite significant differences at the VMC level, DMC has arrived to the same 
value of the energy. Variational parameters differ only slightly for other types of ^He-TJ, 
interaction potentials, except for trimers and tetramers, where the differences are somewhat 
larger. 

We have verified that 1000 walkers are enough for excluding the bias coming from the 
size of the population ensemble used in a simulation. The same conclusion emerged also 
from previous experience in pure and mixed TJ, clusters. The only exception is again the 
'^He2TJ, cluster, for which a larger number of walkers has been necessary: 2000 for the DWW 
potential and 10000 for the TWW potential. 
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In order to eliminate bias coming from the time-step value used in simulations, all cal- 
culations have been performed with several time-steps (At) which assume values within the 
interval 3 x 10"^ — 5 x 10"'^ K~^, where minimum and maximum values in that interval have 
been adjusted depending on the size of the cluster. The ground-state energies for different 
time-steps have then been extrapolated to At — 0. In accordance with the DMC method 
used in this work, which is accurate to second order in the time step,— the extrapolation is 
made with a quadratic function. 

The potential energies and distribution functions have been obtained using pure estima- 
tors.— The asymptotic block size for each cluster and type of ^He-TJ, interaction potential 
has been determined from the behavior of the potential energy versus the block size. For 
most of the clusters, 2500 steps per block have been enough to arrive to the asymptotic value 
of the potential energy. In the case of smaller clusters, having less than five atoms, the block 
size grows to 10000-20000 steps reflecting the low quality of the trial wave functions. Even 
50000 steps in the block have been needed for the TWW potential. However, for blocks of 
this size the same values of pure estimators for different types of trial wave functions have 
been obtained and therefore DMC estimations of these magnitudes are unbiased. 



III. RESULTS 

Ground-state energies of all the studied clusters for two types of ^He-TJ, interaction 
potentials are presented in Table [B Energy in absolute value as a function of the number of 
^He and TJ, atoms is presented in Fig. [2] for the DWW interaction potential. In addition, 
calculations which have been made for several selected clusters with other types of potentials 
are presented in Table [TTl 

As shown in Table [H it takes at least three atoms to form a mixed bound state. The 
weakest bound cluster is ^He2TJ,. However, the value of its binding energy depends strongly 
on the type of the interaction potential. For the DWW, MF and MFmod potentials, it is 
lower than the energy of the ^He2 dimer, which is -1.66 mK, as calculated in Ref. |47l For 



the two weaker potentials (TWW, TY) the obtained energy is within the errorbar equal to 
the energy of the ^He2 dimer. Thus, in the case of the latter two potentials it is unclear 
whether this cluster is stable with respect to the separation in the ^He2 dimer and a free 
T4- atom. On the other hand, for the potentials in which ^He2TJ, is surely bound the 

8 



TABLE I. Ground-state energies (in niK) of the studied ^He-TJ, clusters. In each cell top result 
and bottom result are obtained using the DWW~ and MFmod^ ^He-TJ, interactions, respectively. 



\ Ti 
^He\ 


1 


2 


3 


4 


1 






-56(3) 


-285(5) 






-92(5) 


-355(5) 


2 


-3.0(2) 


-80(6) 


-284(6) 


-635(6) 


-16.6(8) 


-135(9) 


-406(8) 


-819(9) 


3 


-235(10) 


-444(5) 


-790(10) 


-1248(10) 


-290(16) 


-577(5) 


-1007(6) 


-1552(6) 


4 


-790(5) 


-1120(6) 


-1577(5) 


-2145(11) 


-883(4) 


-1318(4) 


-1885(4) 


-2560(6) 


5 


-1664(5) 


-2105(6) 


-2640(6) 


-3310(14) 


-1784(5) 


-2356(8) 


-3014(7) 


-3825(8) 


6 


-2785(7) 


-3335(9) 


-3975(6) 


-4715(13) 


-2925(8) 


-3636(8) 


-4433(10) 


-5320(8) 


7 


-4145(11) 


-4785(7) 


-5510(7) 


-6330(13) 


-4302(6) 


-5121(8) 


-6008(10) 


-7015(8) 




-5725(8) 


-6455(11) 


-7250(10) 


-8155(6) 


8 


-5908(8) 


-6834(12) 


-7818(10) 


-8916(14) 



value of the binding energy varies from -3.0(2) mK to -16.6(8) mK. Such small energies 
result from a huge cancellation between kinetic and potential energies. For instance, in case 
of the DWW interaction potential the potential energy of the ^He2TJ, is Ep = —300(20) 
mK, as determined using pure estimators. The kinetic energy, determined as a difference 
between the total and potential energy is then 297(20) mK. A similar system, ^He2H has 
recently been investigated using an interaction potential which is in between the DWW and 
MFmocl^ and no bound states have been found.— '^'^ This is very much in agreement with 
the present result, since the three times lower mass of hydrogen increases the kinetic energy 



of the system, which already for tritium almost cancels the potential energy. In Ref. 
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other isotopic combinations have been studied as well using hyperspherical coordinates in 
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FIG. 2. The absolute values of the ground-state energies are presented as a function of the number 
of ^He and TJ, atoms, when using the DWW— potential to model the ^He-TJ, interaction. 

TABLE II. Ground-state energies (in mK) of investigated "^He-T^, clusters. 



Cluster 


^He-TJ, interaction potential 


^He 




TWW 


TY 


MF 


2 


1 


-1.3(3) 


-1.4(4) 


-13.9(8) 


3 


2 


-254(12) 


-288(15) 


-522(4) 


3 


4 


-781(11) 


-877(7) 


-1434(8) 


4 


2 


-834(6) 


-886(6) 


-1243(5) 


4 


4 


-1503(10) 


-1633(7) 


-2402(10) 


6 


1 


-2579(6) 


-2605(7) 


-2863(9) 


6 


4 


-3826(10) 


-3990(10) 


-5075(11) 



the adiabatic approximation. One bound state has been found for ^He2T system, with the 
energy of -6.8 mK. Since the authors do not use the same form of the ^He-^He interaction 
potential it is not possible to directly compare this value with our results; however we obtain 
qualitative agreement. 

An extensive search for a bound state of the other mixed trimer ^He(TJ,)2, a possible 
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Borromean state, has been performed. However, for all of the studied potentials and differ- 
ent types of trial wave functions, the energies in the DMC calculations remained positive. 
Moreover, in the course of the simulation, particles were going more and more away from 
each other. This case can be compared with the ^He^He2 cluster which has basically the 
same mass, but interacts with a stronger potential and is likewise not bound.— 

All the other studied mixed clusters with up to four TJ, and eight ^He atoms are bound. 
In absolute value, the energy grows with the addition of both ^He and TJ, atoms, but it 
takes about two atoms to achieve the same increase in binding as with the addition 
of one ^He atom, as can be seen in Fig. [2J Significant cancellation between kinetic and 
potential energies persists for almost all clusters, with total energies being 1% to 20% of 
the potential energy for the largest cluster. The differences in the binding energies for 
different potentials are lowered with the increase of the cluster size and in particular with 
combined lowering of the fraction of TJ, atoms in it. This behavior can be more easily 
seen in Fig. [3] which shows the difference between ground-state energies in calculations 
with MFmod and DWW potentials, divided by the ground-state energy with the MFmod 
potential {EMFmod — Eoas) / EMPmod, as a function of the number of ^He and T| atoms in the 
cluster. Qualitatively, the stability of the clusters, except for ^He2TJ,, is not sensitive to the 
forms of the interaction potentials. However, the relationship of binding strengths between 
different clusters is not always the same for different potentials. For example, ^He7(TJ,)2 
is more strongly bound than ^He6(TJ,)4 for the DWW potential, while for the MFmod one 
this relationship is reversed. Some calculations have been performed using another form of 
the ^He-^He interaction potential^i. The resulting energies are basically equal within the 
errorbars. For example, for ^He2T4, with the HFD-B(He) potential and DWW one obtains 
-3.0(2) mK, and with SAPT— and DWW -3.3(2) mK. The ground-state energy of the larger 
cluster ^HegT; is -5770(10) mK with SAPT, and -5725(8) mK with HFD-B(He). 

The smaller clusters are especially interesting. For example, ^He(TJ,)3^4 has no bound 
two-particle subsystem, while out of all six two-particle interactions in ^He2(T4,)2 only the 



d give a bound subsystem. The latter cluster is similar to 



one between two ^He atoms wou^ 

^He2^He2, which authors of Ref. l50| have called pseudo-Borromean state because of the fact 
that it has also only one bound subsystem. 

The comparison between '^He„(TJ,)m and ^He„^Hem is interesting because both types of 
clusters have essentially the same mass, but differ in the interaction potentials and the statis- 
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FIG. 3. The difference between ground-state energies in calculations with the MFmod^i and 
DWW~ potentials, divided by the ground-state energy with the MFmod potential, IS.E/E = 
{EMFmod — E£)as) / E MFmod-, ^ & function of the number of ^He and T| atoms in the cluster. 

tics. The interaction potential between two helium atoms is the same for all combinations 
of isotopes and is stronger than both the ^He-TJ, and T4,-T4, one. Thus, as expected, the 
clusters "^Hen^Hem are more strongly bound than '^He,i(TJ,)m for m=l,2. However, for more 
than two '^He atoms the fermion statistics comes into play. Thus for example ^He^Hes 4 and 
^He2^He3^4 are unbound, but the clusters ^He(T4,)3,4 and ^He2(TJ,)3^4 have bound states. For 
larger clusters, despite the fermionic nature of ^He, '^He^'^Hem show tendency to be more 
strongly bound than ^He„(TJ,)m, for increasing n, as can be seen in Fig. H] for the DWW 
potential, but is not yet clear for the MFmod. The precise m for which this transition occurs 
can not be predicted because the energies of ^He„^Hem are not calculated for all n and m 
with the HFD-B(He) interaction potential. 

In addition to energy, the structure has also been determined by calculating with pure 
estimators the distribution of particles to the centre of mass of the system p(r) and the 
distribution of the separations between particles P{r). These functions have been normalized 
as j p{r)d^r = N, where N is the number of particles of particular species in the cluster 
and J P{r)dr = 1. Both functions, for the smallest cluster '^He2TJ, are shown in Fig [51 The 
larger plot shows huge differences between distributions of particle separations for different 
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FIG. 4. Comparison of the ground-state energies of ^HeAr(T4,) m clusters for two ^He-TJ, interaction 
potentials with the ground-state energies of ^HeA? '^HeM clusters for different N. The larger figure 



shows results for M = 3, while results for M = 4 are shown in the inset. Axis labe^ 
on both figures. Ground-state energies of ^Hcat '^Hcm clusters are taken from Ref. 



s are the same 



47 and 



51 



potentials. As it can be expected, for all the considered potentials two "^He atoms are on 
average closer than one ^He to a TJ, atom. The average distance between ^He and TJ, atoms 
IS 21.8(5) A in case of the MFmod potential, 47(2) A for the DWW potential and even 
206(15) A for the TWW interaction potential. At the same time < r > between ^He atoms 
is 74(9)A for the TWW potential. This is even larger than the separation of ^He atoms in a 
^He dimer, which for the selected interaction potential is around 50 A. On the other hand, 
in the case of the MFmod potential < r >= 17.0(4) A for the ^He-^He pair (30.4(5) A for 
the DWW potential) which confirms the stability of the cluster. The same floppy behavior 
can be noticed from the distribution of particle distances to the centre of mass, where in 
the case of the TWW potential the probability of finding a particle close to the centre of 
mass is so low that it is not visible in the figure. Correspondingly, the average separation 
of ^He to the centre of mass is < >= 72(6) A and < >= 150(11) A for T|, in case 
of TWW potential, which indicates instability of the cluster. However, in the course of the 
simulation there is no indication that particles are increasing their separation. On the other 
hand, for the DWW/MFmod potential < r^m > is 20.1(5) A /10.2(3) A for ^He and 32(1) 
A/14. 3(4) A for TJ,. Thus, both energetically and from the structure analysis ^He2TJ, is 
predicted to be one of the most weakly bound three-atomic clusters. 

The differences between the distribution functions for the DWW and MFmod interaction 
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FIG. 5. Distribution of interparticle distances P{r) for three different interaction potentials, where 
logscale is used in x-axis. Thinner hnes correspond to the distribution of ^He-^He distances, and 
thicker hnes to ^He-TJ, separations. The figure in the inset shows the distributions of the "^He 
(thinner hnes) and TJ, atoms (thicker hnes) to the centre of mass p(r) for DWW— and MFmod^ 
potential. The distribution functions for the TWW— potential are not visible on the picture. 

potentials decrease with the increase of the cluster size, but are still noticeable in cases where 
the number of TJ, is about the same or larger than the number of the ^He atoms. This can 
be seen in Fig. 6 which presents density distributions with respect to the centre of mass and 
the ^He-^'He and ^He-TJ, separations distribution for the '^He2(TJ,)4 cluster. 

Other small mixed clusters, having in particular 4 to 5 atoms are also very floppy, with 
typical separations between particles of the order of 10 A, for DWW or MFmod potential. 
Figure [7] shows the particle separation distribution for three small mixed four-particle clus- 
ters. As expected, two ^He atoms have on average the smallest separations. The narrowest 
distributions correspond to the most strongly bound cluster ^HesTJ,. With the exchange of 
^He atoms with TJ, ones, clusters become more floppy. It is interesting to notice that in 
the case of ^He2(TJ,)2 and ^He(T4)3 clusters the relationship between particle separations 
distributions for ^He-T| and T^-TJ, pairs is the opposite. 

When TJ, is added to the helium cluster ^HeAr, as a consequence of its smaller mass it 
tends to stay further away from the centre of mass of the system than ^He atoms, as can 
be seen in Fig. [HI This is completely analogous to the behavior of '^He in ^HeAr^He system. 
The only exception being that TJ, is on average further away from the centre of mass in the 
"^HejvTJ, cluster than '^He atom in the "^HeTv^He cluster, which is a consequence of the weaker 
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FIG. 6. Distribution of of the ^He (thinner hnes) and TJ, atoms (thicker hnes) to the centre of 
mass p{r) for DWW— and MFmod^ potentials. The interparticle distance distribution P{r) is 
shown on the figure in inset. Thinner hnes correspond to the distribution of ^He-^He distances, 
and thicker hnes to ^He-TJ, separations. P{r) for T4,-T4, are very close to P{r) for ^He-TJ, and are 
thus not shown for better clarity. 



l0^P(r)(A:') 




FIG. 7. The distribution of the interparticle distances for three mixed tetramers. 



^He-TJ, interaction potential. This can be seen by comparison with the density distribution 
of ^He in Ref. 



52l . For = 3 there is still appreciable probability for TJ, to be close to the 
centre of mass, while for = 8 this probability is already very low, with < r^m >= 7.33(1) A 
and root-mean-square deviation Arcm = a/ < > — < fcm >^ of 2.2 A. For larger clusters, 
it is expected that TJ, will be completely pushed to the surface forming the so-called Andreev 
states,— like ^He in mixed ^He-'^He clusters.— 
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FIG. 8. Density profiles of tlie '^HeTvTJ, system, for N=3,5,8. Thin and thick lines stand for ^He 
and TJ, atoms, respectively. Left scale is for ^He and the right one for TJ, density distributions. 



Having determined the average squared interparticle separations in DMC by pure esti- 
mators, it is interesting to investigate if any of the smallest clusters would fit the definition 
of the quantum halo state. For three-body systems a useful definition has been given in Ref. 
271 . Quantum halos should obey the condition (p^) /p^ > 2, where p is the average radial 
coordinate and po a scaling parameter, defined as 



mp 



M ^ 



miUikiri - Vk 



mpl 



1 v-^ 



(7) 



i<k i<k 

In the previous equation, m is the arbitrary mass unit, rrii is the mass of atom z, M is the 
total mass and Rik is the two-body scaling length of system i and k. For two-body systems 
Rik can be easily determined as a classical turning point, from E = V{r). In order to make a 
rough qualitative prediction we here assume that on average each pair contributes with the 
same amount to the total binding energy. The same assumption has been made in Ref. I29 . 
where clusters of particles of the same type have been considered. Here, this assumption is 
obviously not fully correct, but we have verified that the final conclusions whether a cluster 
is a candidate for halo state or not, remain the same if we assume absolute average potential 
energies of some pairs to be larger and others smaller. Finally, since the classical region is 
defined by having positive kinetic energy, we estimate the classical radius by equaling the 
total energy divided by the number of pairs with the corresponding potential energy of each 
pair. Having thus found all the Rik we calculate the scaling parameter po from Eq. ([?]). 
The only bound mixed trimer '^He2(T4) fits the definition of the quantum halo state since 
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(p^) /Po = 13.6. The result is insensitive to different types of ^He-TJ, interaction potentials. 



The authors of Ref. 



29| have used the same approach to estimate the halo condition for 



clusters of four and five particles. If we extend Eq. ([7]) to larger clusters we find that 
both ^He2(T4)2 and ^He(T4)3 can be considered quantum halo states in case of the DWW 
■^He-TJ, interaction potential, since (p^) /pi is 2.32 and 2.04, respectively. This approach is 
qualitative, but it confirms the extremely weak structure of these systems and motivates 
their further investigation. 



IV. CONCLUSIONS 

The ground state of small mixed clusters composed of helium and spin-polarized tri- 
tium has been investigated by quantum Monte Carlo simulations. Significant differences in 
binding energies are obtained for different ^He-TJ, interaction potentials, especially for the 
smallest clusters. Nevertheless, the conclusions concerning stability limits are insensitive 
to these differences, except for the trimer. The mixed dimer does not exist, and the only 
bound mixed trimer is the ^He2T^ one. The latter result is obtained for several ^He-TJ, 
interaction potentials. However, in the case of the weakest potential (TWW) the cluster 
^He2TJ, is at threshold of separating into a dimer and a free particle and its stability could 
not be determined with certainty. Due to its large size, this cluster has most of the proba- 
bility outside the classically allowed regions of space, classifying it as a quantum halo state. 
Another possible mixed trimer ^He(TJ,)2, a candidate for the Borromean state is not bound. 
Our results further indicate that two mixed tetramers ^He2(TJ,)2 and "^He(T4,)3, could also 
be considered quantum halos. The largest clusters here analyzed are also very weakly bound 
and floppy. 
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